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Abstract 

Scenarios where gravitinos with GeV masses make up dark matter are known to be 
in tension with high reheating temperatures, as required by e.g. thermal leptogenesis. 
This tension comes from the longevity of the NLSPs, which can destroy the successful 
predictions of the standard primordial nucleosynthesis. However, a small violation of 
matter parity can open new decay channels for the NLSP, avoiding the BBN prob- 
lems, while being compatible with experimental cosmic-ray constraints. In this paper, 
we propose a model where matter parity, which we assume to be embedded in the 
U(1)b-l g au g e symmetry, is broken dynamically in a hidden sector at low scales. 
This can naturally explain the smallness of the matter parity breaking in the visi- 
ble sector. We discuss the dynamics of the corresponding pseudo Nambu-Goldstone 
modes of B — L breaking in the hidden sector, and we comment on typical cosmic-ray 
and collider signatures in our model. 



1 Introduction 



Rapid proton decay in supersymmetric extensions of the standard model is usually avoided 
by the introduction of R- (or matter-) parity. However, it turns out that in scenarios with 
gravitino dark matter with masses in the range 10 to a few 100 GeV, a small violation 
of this symmetry can be desirable since this would reconcile thermal leptogenesis as a 
source for the baryon asymmetry in the universe with standard primordial baryogenesis [lj 
(for a recent exhaustive study of small bilinear i?-parity violation see |2]). The tight 
upper bounds on the corresponding parity-violating couplings, stemming from laboratory 
experiments OH] and the non-observation of gamma-ray lines from gravitino decays at 
Fermi LAT [5J, then raise the question how the presence of such small couplings can be 
understood theoretically. In this paper we consider the possibility that the small violation 
of matter parity, which we assume to be embedded into the U(l) B _ i gauge symmetry, is 
linked to the condensation of strongly interacting hidden sector fields at low energies. 

Thermal leptogenesis [6j via the out-of-equilibrium decay of heavy right-handed neutri- 
nos is one of the most promising models for the generation of the baryon asymmetry in the 
universe. Furthermore, beside explaining the observed excess of baryons over anti-baryons, 
the introduction of right-handed neutrinos can also naturally explain the smallness of the 
neutrino masses; this is known as the see-saw mechanism 0H]. Leptogenesis as well as the 
see-saw mechanism favor very high masses for the right-handed neutrinos, and it is well 
known that especially thermal leptogenesis requires that the mass of the lightest right- 
handed neutrino should be larger than 10 9 GeV which gives also a strong lower limit 
on the reheating temperature. 

In supersymmetric extensions of the standard model, gravitinos play a distinct role due 
to their Planck-suppressed couplings to the observable world, and they are an excellent 
candidate for being the dark matter observed in the universe. It is well known that, 
given a high reheating temperature Tr ~ 0(1O 9 -1O 10 GeV) and gravitino masses around 
rri3/2 ~ 0(10-100 GeV), the thermal relic density of gravitinos reproduces the correct dark 
matter abundance |10| . in good agreement with thermal leptogenesis. 

However, the above scenario is not free of problems: The next-to-lightest supersym- 
metric particle (NLSP), which is often the stau, typically decays into the gravitino with 
lifetimes that are of the order of minutes or days. This is generically in conflict with the 
successful predictions of standard primordial nucleosynthesis (BBN), since the NLSP decay 
releases large amounts of electromagnetic or hadronic energy into the primordial plasma, 
destroying some of the fragile elements and changing their abundance [TT|I12|. Further- 
more, if the NLSP is charged, it could form bound states with 4 He and 8 Be, leading to a 
catalytic over-production of 6 Li and 9 Be [T3],ITI]. 

Different solutions to the above problem were proposed: If in some special scenarios 
the sneutrino or the stop is the NLSP, their late decay would not affect the predictions 
of BBN strongly [T5jlT6]. In case of stau NLSPs, a large left-right mixing could lead to 
suppressed relic abundances [IT]. Entropy production between thermal freeze-out of the 
NLSP and BBN could dilute the NLSP and gravitino abundance enough to evade the 
BBN constraints [T8}IT9]. The NLSP mass could be nearly degenerate with the gravitino 



2 



mass, which reduces the released energy [20]. Or the NLSP could decay into additional 
hidden sector states before conflicting with BBN |2~T||2"2"]. Lastly, as mentioned above, in 
cases where i?-parity is weakly violated, the NLSP can decay into standard model particles 
before the onset of BBN pQ. 

Theoretically, the existence of i?-parity is not automatic. If it is a global symmetry it is 
expected to suffer from quantum gravitational 'wormhole' effects [23J24J. Therefore, it is 
necessary to assume an embedding into a local symmetry, which is spontaneously broken at 
a high scale. However, the requirement of anomaly cancellation strongly disfavors gauged 
XJ(1) R proposals. Another possibility is to understand discrete i?-symmetry [25j in the 
framework of higher- dimensional constructions. The gauging then corresponds to local 
coordinate rotations in the extra dimensions, and the i?-parity of the effective theory is a 
relic of the higher-dimensional Lorentz group. In string theory such discrete i?-symmetries 
may appear in specific compactifications [26], however, not generically. 

Alternatively, proton stability in supersymmetric setups can be guaranteed by matter 
parity. This is natural in models which incorporate the see-saw mechanism, since the 
required Majorana mass Ms of the right-handed neutrinos violates B — L by two units, 
thereby breaking U(1) B _ L to its matter parity subgroup. Note that U(1) B _ L is the only 
anomaly-free Abelian extension of the standard model, and that it is canonically included 
in SO(10) GUT proposals. Therefore, we consider the 'matter parity' to be the Z 2 subgroup 
of the anomaly-free U(1) B _ L gauge symmetry throughout this paper. 

As argued above, the violation of i?-parity or matter parity is bound to be very small. 
This can be achieved by e.g. connecting i?-parity breaking with B — L breaking, via non- 
renormalizable operators [1] or at loop level in the context of string theory [27J. 

In this paper, we discuss a scenario where the small matter-parity breaking is related 
to the condensation scale A of an asymptotically free hidden sector gauge group factor. 
This scale can be naturally much lower than the U(l) B i breaking scale Ms- If the hidden 
sector is charged under B — L, the condensate breaks matter parity at low scale. It then 
dominantely couples to the right-handed neutrinos, thereby inducing a sneutrino vacuum 
expectation value and consequently small bilinear matter parity violation [28] of the order 
A 2 /M s , which can lie in the phenomenologically preferred regime. 

The paper is organized as follows: In the next section we will introduce the hidden 
sector model and discuss its coupling to the visible sector; in section 3 we will discuss 
phenomenological aspects of the model, namely gravitino LSP and stau NLSP decay, ther- 
mal production of hidden sector particles and their decay; section 4 is devoted to further 
discussions and in section 5 we will conclude. 

2 A Model for Dynamical Matter- Parity Breaking 

We consider an extension of the supersymmetric standard model with gauge group 

G = SU(3) C x SU(2) L x U(l) y x V{1) B _ L x SU(2) hid , (2.1) 
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where U(l) B _ L will be broken spontaneously to matter parity at the high scale Ms, and 
SU(2) hid is the gauge group in the hidden sector, responsible for the low-energy matter- 
parity breaking. 

We assume that light neutrino masses are generated in the standard way by the see-saw 
mechanism after the breaking of \J(1) B _ L at a high scale M5. Large Majorana masses for 
the right-handed neutrinos N c are generated by a singlet field S in our model, 

Wsee-saw =h\f + ^\?SN?N? , (2.2) 

where we have taken the basis where the Majorana mass matrix for N c is diagonal. We 
assume (S) = Ms- Then, the N? have Majorana masses, Mj = Af x Ms for i — 1 — 3, 
respectively, and we define M\ < M2 < M3 (we will assume that M 3 ~ Ms below). 
For successful thermal leptogenesis, we require Mi > 10 9 GeV. Since U(l) B _ L is gauged 
initially, the corresponding Nambu-Goldstone (NG) mode will be absorbed by the gauge 
field, which then decouples from the low-energy effective theory. The field S necessarily 
has B — L charge —2, hence a global discrete 7L 2 subgroup remains unbroken during this 
process^, the 'matter parity'. We also assume the presence of a field S with B — L charge 
+2, required for anomaly cancellation, without discussing the dynamics that give rise to 
the condensation of S and S in this paper. 

The remaining matter-parity is then broken dynamically by the condensation of quarks 
in the gauge group SU(2) hid at the low scale A <C Ms- The hidden sector contains two 
doublet quarks Q", Q 2 with B—L charge +1/2, and two doublet quarks Q3, Q2 with charge 
—1/2. Furthermore, we require the existence of five neutral singlets Zi 3 , Z14, Z 23 , Z24 and 
X. Here and in the following a, (3 — 1,2 are SU(2) hid indices. The low-energy degrees of 
freedom are the antisymmetric combinations [29J 

fy = -fy = ~Q?e ia , (2.3) 

with convention Q ia = e a /3<2f , where e a p is totally antisymmetric and ei 2 = 1. 
We consider the effective superpotential (cf. [291130]) 

W dyn = X (Pf (yd - A 2 ) + A (Z 13 V 13 + Z U V U + Z 23 V 23 + Z 2 ,V 2i ) , (2.4) 

where Pf (Vij) = V12V34 + V14V23 — V13V24 is the Pfaffian of the antisymmetric matrix Vij and 
the X is some composite state. The accidental symmetries of the model will be discussed 
in Section HI 

The superpotential (12. 4|) possesses a unique supersymmetric minimum (which is in con- 
trast to [30]), with non- vanishing vacuum expectation values for the two charged effective 
mesons, 

(V ) - (V ) - A ^ = ^ = ^ = ^ = ° ' (2 <n 



^^The generator of matter parity combines B — L with hypercharge, P = {- 



\B-L+4Y 
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Here, we assumed that (V12) = (V34) is a consequence of the soft mass terms for the Vy, 
induced by SUSY breaking. Since Vi 2 and V34 have B — L charge +1 and —1, respectively, 
we conclude that matter parity is broken dynamically at the low scale A in the hidden 
sector. 

The only unsuppressed and renormalizable interaction in the superpotential, connecting 
hidden and visible sector, is given by the term 

W D -fiQsQ^Ni . (2.6) 

Here, fi with i = 1,2,3 are free parameters that we will fix below. In the vacuum ( 12. 51) 
this becomes a linear term for the right-handed neutrino multiplets, which together with 
the mass term in Eq. (12.21) implies a non-vanishing vacuum expectation value for the 
corresponding sneutrinos, 

This vev generates matter parity violating couplings between the standard model fields in 
the superpotential [28J , which are of the order 



W DO (—} LH U + O ( — ^ LLE C 

\M S ) \M S M P ) 

+ (tt^yt-) l Q dc + (ttKt-) D c D c U c , (2.8) 
\M S M P J \M S M P J v ; 

coming from gauge invariant terms like Mp 1 N c LLE c in the superpotential. Thus the 
matter-parity breaking is mostly bilinear in our model, originating from the Yukawa cou- 
plings in (12. 2p . Its scale is related to the condensation scale of the hidden sector gauge 
group. The effects on the matter-parity violation depend on many unknown Yukawa cou- 
plings and we will assume a very simple situation where fa < 1 and fx = f 2 — in this 
paper, to show the presence of a consistent parameter region in the model. Note that we 
consider that the N° decay is relevant for leptogenesis. Since fi = 0, the N° decays into 
the hidden sector are suppressed and hence leptogenesis is not affected by the presence 
of the hidden sector. Furthermore, the matter-parity violating couplings in Eq. (12 .8p are 
small enough to not wash-out the lepton asymmetry again j31] for the parameters we are 
interested in pQ. 

When the interactions (12.61) between hidden and visible sector are neglected, we have a 
global XJ(l) B L symmetry in the hidden sector, since the corresponding gauge-symmetry 
is already broken by the condensation of S and S. This global U(1) B _ L symmetry is 
then broken by Q a Q a condensations at the scale A, producing a massless pseudo NG 
multiplet in the hidden sector. Through the interactions (12. 6p . this multiplet would receive 
a supersymmetric mass m^Q = f£j^ ~ C(keV)Jl In addition, there are hidden sector 



2 The true NG boson of the B — L breaking is a linear combination of this pseudo NG boson and the NG 
boson arises from the high-energy B — L breaking. This true NG boson is nothing but the one absorbed 
in the B — L gauge field. 
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multiplets with a mass ~ A. However, after supersymmetry breaking, the soft mass terms 
raise the masses of the modes in the pseudo NG multiplet. For the pseudo NG boson this 
yields a mass (cf. app. [B]) 

m „ „ f , ^C^g K 1M MeV/ 3 ( J=^) 1 ( I5 A^ 7 ) ~* , (2.9) 

whereas the tree-level masses of the fermion partner ip and the radial scalar pseudo NG 
component p are given by 

( m ln\ 

m^p ~ %/2 + ^ I ) ' m p ^4m 3/2 , (2.10) 

where we assumed that SUSY breaking is mediated to the hidden sector only via gravity 
mediation. Below we will always assume that > m 3 / 2 , since we are interested in 
gravitino dark matter. 



3 Phenomenology 

In this section we will discuss the decay of the gravitino and the NLSP, which are general 
features induced by the matter-parity violation, as well as the dynamics of the pseudo NG 
modes of B — L breaking in the hidden sector, which are special to the model we proposed. 
The resulting bounds on the parameter space are summarized in Fig. [1] for two reference 
scenarios. 

3.1 Gravitino and NLSP decay 

As stated above, the matter-parity violating effects on the MSSM sector are mainly induced 
by the vev of Nf and hence dominated by the bilinear term 

W D lnUE u , (3.1) 

which is generated in the superpotential after condensation of the hidden SU(2) hid gauge 
group as described in the previous section. Here, following Eqs. (12. 2\ and ( 12. 7p . the size of 
the matter-parity violation is given by 

* = f ^W 3 ^ w " GeVf > (#) (lo4v) 2 (to)" • (3 ' 2 > 

where we assumed that / 3 < 1 and f\ = f% = 0, as mentioned above. 

For the phenomenological analysis, it is convenient to perform a rotation between the 
higgs- and lepton-doublet superfields: 

H d — >■ H d — EiLi , Li — > Li + EiH d , (3.3) 
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where £j = fii/fi, and fi denotes the MSSM //-term. The rotation eliminates the bi- 
linear term in Eq. (13. ip . but it generates the matter-parity violating trilinear terms 
\Xij k LiLjEl and X'^LiQjDl with coefficients given by X ijk = and A- jfc = -ejif^, 

respectively. In general, bilinear i?-parity breaking terms such as (13. f I) imply a small 
but no n- vanishing vev for the left-handed sneutrinos z/j. After the rotation (13. 3p it 
reads (z>j) ~ — e^v^ (1 + ^ 2 /rh 2 Li ), where fh 2 Li denotes the soft slepton masses and soft im- 
parity breaking terms were ignored^]. Note that the operator ^\'lj k U^DjDl, which finally 
leads to proton decay, is not generated in this way. Nevertheless, due to the Planck- 
suppressed operator Mp 1 N C U C D C D C , it is generically present in the superpotential with 
coefficients A" ~ 0(A 2 /MpM%), cf. Eq. (12. 8p . Constraints from proton- stability require 
^lifc^iifc ~ 10~ 27 [32| . and a sufficient condition to satisfy them is A < 10 6 GeV, where 
we assumed that M 3 ~ 10 16 GeV and f 3 = 1. This tight constrained will be in general 
relaxed depending on the flavor structure of \'[- k and and on the value of / 3 , and we 
will assume this to be the case below@ 

Most importantly, the above violation of matter parity induces the decay of the grav- 
itino [1]. The corresponding decay widths for the two-body decay into is given by |33| 



3 

1 Tn 

m-m^i») = ^W^, (3.4) 



where U^ u parameterizes the photino-neutrino mixing, and M P = 2.4 x 10 iS GeV denotes 
the Planck-mass. Assuming that SUSY breaking only induces matter-parity conserving 
soft-masses at tree-leve0, it can be approximated as |2] 

s w M w -Mg _,,„_,,/■ M B 



where ( = — £ i an d v d — v cos ft = (H^) denotes the higgs vev. We made use of the 
GUT relation M-^/M^ ~ 2, where Mg and are the bino and wino masses, respectively. 
From this, we find for the inverse decay width 

r,^ . ^ « 4 x ^ . (^)- 2 ( I J£ fe7 ) 2 (^) 2 , (3,) 



3 Vanishing soft bilinear i?-parity breaking terms can in fact only be demanded at one particular scale. 
In the generic case without fine tuning bilinear i?-parity breaking has 9 parameters, whose role for gravitino 
and NLSP decays was analyzed in detail in [2]. Here we ignore phenomenological complications due to 
the presence of soft bilinears, since in this paper we are mainly interested in the theoretical origin of the 
supersymmetric term (|3.1[) . 

4 Note that there is also the Planck-suppressed operator QQQL in the superpotential, which leads to 
proton decay and which is not forbidden by the B — L symmetry. We will assume that this operator is 
suppressed by some unspecified mechanism, and note that it is not significantly regenerated by breaking 
of matter-parity. 

5 Note that the above rotation (|3.3|) regenerates Bh u hd, rh 2 Li l*li and fh 2 d h* d hd, which then determine 
the photino-neutrino mixing angle [2]. 
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with e = a/S £ i — as l° n g as £ i) £ 2 ^ £3- Other decay channels like ip 3 / 2 —> W ± £ T 
and ip 3 /2 — >■ and decay channels with the hidden sector pseudo NG boson in the 

final state, also contribute to the overall decay widths. However, since the branching 
ratios are of comparable size [1] for the different two-body channels into SM particles, we 
will concentrate on the decay into gamma-ray lines. This decay channel is constrained 
by gamma-ray observations, which will give strong bounds on the matter-parity violating 
parameter e 3 in the region of interest. 

The decay of the NLSP, which we take to be the lightest stau T\ for definiteness in the 
rest of the paper, is mediated by the term LiLjE^- For simplicity we will assume that t\ is 
dominantly right-handed, otherwise LiQjD c k would induce also decays into hadronic jets, 
see Ref. pQ. The corresponding decay width is dominated by decay into v^t and v T fi and 
given by 

r(ri ->■ u^r, u T fi) ~ 2 x -^|A 23 3| 2 ^n , (3.7) 

107T 

where we neglected fermion masses, and A2 3 3 = £3^3 — £2^33 ~ — 0.1 (tan /3 / 10) £2- From 
this we obtain for the corresponding decay time 

rfi ~ , x ir' s (^) ' , (3,) 

from which it follows that the decay of the NLSP can easily happen before BBN. Below 
we will always assume that e 2 — £3 — £ for simplicity. 

As mentioned above, a number of astrophysical observations constrain the lifetime of 
the stau and the gravitino, which translates into bounds on the matter-parity breaking 
parameters. Firstly, the t% has to decay before BBN, with a lifetime shorter than 2 x 
10 3 s, to avoid catalytic overproduction of 6 Li [T3"|IT4"]. This implies the lower limit of 
roughly e > 10~ 13 . Furthermore, the gravitino decay into gamma-ray lines is limited 
by observations of the Fermi LAT satellite to lifetimes r^^^ > 10 29 s, yielding an 

upper limit of e < 6 x 10- 9 (m 3/2 /100 GeV)" 3 / 2 in the mass range where the bounds were 
evaluated, 60 GeV < m 3 / 2 < 400 GeV |34] . The corresponding limits on the condensation 
scale A are summarized in Fig. [T] for two reference scenario of parameters. In these plots 
we also show the strong lower bounds on the scale A and upper bounds on f'3 that arise 
from the late decay of the pseudo NG modes of the hidden sector, which we will discuss 
next. 



3.2 Thermal production of hidden sector particles 

We will now discuss the cosmology of the hidden sector containing the condensing SU(2)hid 
gauge group. Even when conservatively assuming that reheating only affects the MSSM 
sector, hidden sector particles are produced in the early universe by scattering of MSSM 
particles. This can potentially lead to overclosure and can reintroduce problems with BBN. 
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Figure 1: Summary of constraints on condensation scale A, as function of the free parameter 
/3, for two reference scenarios with Ms — 10 16 GeV (left panel) and Ms = 10 13 GeV (right 
panel). We assumed m 3 / 2 = 120 GeV, y = 10 _1 and a reheating temperature of Tr — 
10 9 ' 8 GeV, which yields the correct relic density for the thermally produced gravitino [10]. We 
show the bounds coming from catalytic production of 6 Li during BBN by late decaying staus 
(BBN stau J14], for stau masses around 150-200 GeV), from the late decay of the fermionic 
partner of the pseudo NG boson of matter-parity breaking (BBN NG, from Ref. \12) , assuming 
branching ratio into hadrons ~ 1), and from the non-observation of gamma-ray lines from 
gravitino decay by Fermi LAT. For comparison, we show also the region where the pseudo 
NG boson relic density exceeds the observational limit f2/i 2 ~ 0.11 (OP NG). As apparent 
from the plots, bounds related to the pseudo NG modes become more sever for smaller values 
of Ms- We always assumed fx = f 2 = 0. 



From (12. 4p it follows that the hidden sector contains 10 chiral multiplets with a mass 
at the condensation scale, ~ A, as well as one chiral multiplet that remains light and is 
the pseudo NG mode of the breaking of the accidental global B — L symmetry in the 
hidden sector. As mentioned above, the pseudo NG boson a acquires a small mass term, 
suppressed by M s , due to the coupling to the B — L breaking visible sector, see (12. 9p . 
whereas the fermionic partner if) and bosonic partner p of the pseudo NG boson get masses 
around the gravitino mass 777.3/2, see (I2.10p . 

At temperatures above the condensation scale, T ^> A, the production of hidden sector 
particles is dominated by the Ms-suppressed dim-5 operator 

W D -j-r JzQzQiLHu , (3.9) 

which is generated after integrating out N£, and by the M| suppressed dim-6 operators 
coming from integrating out the heavy B — L vector- multiplet, which gives operators like 

T ^ (gigl + Q\Q,2)L ] L 

K D (3.10) 

in the Kahler potential K [33]. Here we have defined y = max^^, , h^), and we will 
ignore the flavor structure for simplicity. Note that the term (I3.9P violates the U(1) B L 
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symmetry explicitly, but since this interaction between hidden and visible sector is never 
in equilibrium, the lepton asymmetry produced by leptogenesis is maintained. 

In the first case, the matrix element for 2 — > 2 scattering processes like lh u — > Q^Qi 
is of the order \Ai\ 2 ~ y 2 f 2 s/M 2 , where s denotes the center-of-mass energy. Following 
Ref. |36| . the yield Yq. = uqJ S of hidden sector quarks Qi, where S is the entropy density, 
is given by 

where we have used g* ~ 915/4 as the number of effective relativistic degrees of freedom of 
the MSSM. Similarly, in the second case, the operator in (13 . 1Q[) induces 2 — > 2 scattering 
processes like /*/ — > Q\Q\, and the corresponding matrix element is of the order \Ai\ 2 ~ 
s 2 /Mg. The corresponding yield is given by 

T 3 R M P T R \ 3 / M s V 4 fo,n\ 

Qi ~^MW 2 ~ U^GevJ iTo^GevJ " (3 " 12) 

The above processes are most efficient at the highest temperature, Tr. The energy 
density produced in the hidden sector can be hence estimated to be pns ~ ^r^hs> where 
^hs — Yq V S denotes the particle density in the hidden sector shortly after reheating. 
Assuming that particles in the hidden sector quickly thermalize, this corresponds to a 
temperature Ths in the hidden sector that is given by (tt 2 / '30) ] T^ S = pns [37] . From 
this it follows that the temperature in the hidden sector shortly after reheating is Ths ~ 

(^vsV^H^' ) ) 1//4 (^ / a V ) 1//4 ^ fi; ' wnere 9bs an d 9ys denote the corresponding effective degrees 
of freedom in the hidden and visible sector at high energies. 

Depending on the details of the symmetries of the hidden sector, some of the heavy 
hidden sector states with masses around the condensation scale ~ A can be stable, which 
could lead to overclosure of the universe. However, despite their large masses, which can 
exceed the unitarity bound A < 10 5 GeV [38], sufficiently strong annihilation of the hidden 
sector states is in principle possible due to the large temperature difference between hidden 
and visible sector, which reduces the relic density by a factor of order Ths/T/s ~ (^s^) 3 ' 
see Ref. [39] • However, a detailed study of this annihilation process is beyond the scope 
of this paper. Subsequently we will assume for simplicity that all these potentially stable 
heavy hidden sector states efficiently annihilate into the pseudo NG modes, and we will 
consider the subsequent effects of these NG modes only. Alternatively, depending on the 
underlying symmetries, the heavy hidden sector particles could directly decay into Standard 
Model particles before BBN, thus avoiding the overclosure problems from the beginning. 
We will discuss this situation shortly in section H] below. 

To estimate the hidden sector temperature and density of pseudo NG modes at 
later times, one can make use of the fact that the entropy densities, which is given by 
'S'hs/vs — (27r 2 /45)(y'HS/vs^HS/vs' ^ n ^ ne hidden and visible sector are conserved sepa- 
rately [37]. >From this it follows that today T H s ~ (5 , vs/5 , Hs) 1 ^ 3 (^ / g I i V ) 1/ ' 4 2vs; where g V s an d 



10 



<7hs denote the corresponding effective degrees of freedom of the hidden and visible sector, 
now at low energies. Noting that the particle density today is given by n^s ~ ^hs/^^hS' 
this gives an effective particle yield of Y^q = nns/'Svs ~ 0.2(Yg. v ) 3 ' 4 . The corresponding 
relic density of the different pseudo NG modes is given by 

nE>» = 2.74 x 10- (J£gS_) O , (3.13) 

and we take for definiteness 4F a uv = 4F p uv = 2Y^ Y = Y$f. 

At temperatures below the condensation scale, 100 GeV < T < A, hidden sector 
particles are predominantly produced via the renormalizable operator 

W D yf 3 ^-V u LH u . (3.14) 

Here, the decay of particles in the thermal bath into the hidden sector particles is relevant, 
like e.g. decays of the higgs h° — > ipv. The corresponding relic density is given by [36] 



W gV % y™5 U00GevHl0 6 GeVy' \10 16 GeV 

(3.15) 

where Tj ~ l/87T(yf 3 A/M 3 ) 2 mj denotes the approximate decay widths of particles j in the 
bath into the hidden sector modes, and we used ~ 100. This can dominate over (13. 13[) 
when A is large. 



3.3 Decay of the hidden sector pseudo NG modes 

Finally, we discuss the decay of the pseudo NG modes. The fermionic pseudo NG boson 
partner ijj mixes with the neutrinos z/j and inherits their SM interactions. This mixing is 
induced by the term (12. 6ft . and given by (for details see app. [A]) 

which is suppressed by the B — L breaking scale Ms- Here, the index i refers to neutrino 
flavor eigenstates. Using ( 13. 2p . the size of the mixing angle can be expressed in term of the 
matter-parity breaking parameter like 

*< = 12 x 10 - 13 (™) (ra)" 1 (io4v) _1 ■ < 317 » 

If > Mz, the two-body decay if) — > Z°vi is kinematically allowed, and the corresponding 
decay width is given by 



^ l) - 16*8^ Ml [ 1+ ml ml 
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— (0.03 sj [ w -u) [ lbQGeY ) M 2\^+ ^ Z ml){ 1 ml)- 

(3.18) 

For masses around m^ ~ 150 GeV and mixing angles 9i ~ 10~ n this gives inverse decay 
widths of the order of r ~ 0.03 s, which is well before and save for BBN |12j . However, 
when 9i is smaller, BBN bounds become very relevant |§ The lifetime of i/j is further 
reduced by the decay channels ijj — > W whose decay width is obtained from ( 13. 1 8[) by 
r^p^T = ^c^yT^zo^. and substituting M z —> M w ±. Note, that if m^ < M w ±, the 
much slower three-body decays like ip — > e + e~z/ would become important. They would 
give rise to much stronger lower bounds on the mixing angle 9 than in the regime where 
two-body decay is allowed. 

The bosonic pseudo NG boson partner p can perform the two-body decays like p — > h u v, 
if kinematically allowed, as follows from the effective operator (I3.9p . This requires that 
M~ h < m p . The lifetime of p is typically much smaller than the lifetime of ip, since its 
decay is not further suppressed by gauge couplings. When the above decay is kinematically 
forbidden, three-body decays like p — > bbv become relevant. However, this depends on the 
details of the MSSM spectrum, and we will conservatively assume that the decay of p is 
always faster than the decay of ijj. 

The pseudo NG boson a in principle also decays into SM particles. The interactions 
with the visible sector are at low energies mainly mediated by the operators V^LiH u , 
and therefore the interaction is suppressed only by a factor A/M3. However, since 
the particle interacts with the visible sector only via a matter parity violating operator, 
final states for its decay can only be odd under matter parity. But no matter parity odd 
final states with integer spin number exists below the mass of the lightest gaugino, and 
hence the particle is essentially stable and one has to require that f2jj R,uv /i 2 < 0.1 to avoid 
over closure 

Lastly, we note that the invisible decay of gravitinos into the pseudo NG boson and neu- 
trinos, which is kinematically allowed, happens with lifetimes around r ~ 10 32 s (10~ l2 /9) 2 
for ~ 100 GeV masses, and can hence be neglected. Furthermore, the decay f — > Tip, which 
is also in principle possible, is double suppressed by the small matter-parity breaking p/p, 
and hence negligible with respect to the decay mode f — > tv. 

The bounds on the parameter space of our model are summarized in Fig. [T] for two 
reference scenarios. We show bounds from the non-observation of gamma-ray lines with 
Fermi LAT, where we have assumed that gravitinos make essentially all of dark matter 
today. Furthermore, we show bounds from BBN that stem from the late decay of the stau 
via the induced R-parity violation, and bounds from BBN that come from the late decay 

6 Note that the decay into gravitinos, ip — > aV'3/2, is typically very slow with lifetimes of the order of 
r > 10 8 s or larger, and can be safely neglected. 

7 In fact, the decay rate is doubly suppressed by the small matter-parity violation, and loop suppressed 
since the decay involves the higgs-multiplets, e.g. T(a —> 77) ~ a 2 (yfsA/M 3 fl/fj.) 2 m^/fi 4 ~ O(10 42 s) -1 . 
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of the fermionic and bosonic pseudo NG boson partner. Lastly, we also show the region 
where the relic density of the pseudo NG boson would exceed the observational limit. 

4 Discussion 

As shown in Fig. HJ cosmologically viable models with high reheating temperatures and 
gravitino dark matter can be found for large enough XJ(l) B L breaking scales Ms > 
£>(10 13 GeV) and SU(2) hid condensation scales A ~ O(10 3 -10 8 GeV). We also assumed 
suppressed couplings between hidden sector and the lighter two right-handed neutrinos, 
adopting j\ = f 2 = for simplicity. In these models, the communication between the 
matter-parity breaking SXJ(2) hid sector and the visible sector is suppressed by the mass of 
the heaviest right-handed neutrino, M3 ~ Ms- As an important consequence, the mass 
of the pseudo NG boson of the breaking of the accidental global B — L symmetry in the 
hidden sector, as well as its relic number density, are small enough to yield a relic den- 
sity in agreement with the observations, so that the pseudo NG boson constitutes only a 
subdominant component of dark matter. Furthermore, the relic density and the lifetime 
of the pseudo NG boson partners can be small enough to be in agreement with the tight 
BBN constraints. Note that requiring that the lighter two right-handed neutrinos do not 
change the above phenomenology implies upper limits on fi and f 2 that are at least of the 
order fi < f^Mi/M^, since the interactions between hidden and visible sector depends on 
the ratio fi/Mi. 

Due to the strong BBN bounds on the relic abundance and the lifetime of the pseudo 
NG modes of B — L breaking, large values of matter-parity breaking are favored in our 
model. This implies both a lifetime of the gravitino closer to the experimental limit, as well 
as shorter decay lengths of the NLSPs. For values of matter-parity breaking close to the 
Fermi LAT limit, we estimate the decay length of the stau to be of the order of (9(100 m). 
It is very interesting, that this is in principle accessible by the LHC |3D]. However, the 
quantitative details depend e.g. on the flavor structure of the stau couplings, and we refer 
to Ref. [2j for a detailed discussion. 

In general, the low-scale breaking of discrete symmetries can induce serious cosmological 
domain wall problems, which arise after the completion of inflation and thus are never 
diluted pi] . In the presented scenario these problems can be avoided in two different 
ways. 

Firstly, estimating the maximum temperature Ths of the hidden sector immediately 
after reheating, by equating pns ~ ^vs^vs^ ~ ^~hs> snows that the hidden sector temper- 
ature after inflation does never exceed the breaking scale A in a large part of the allowed 
parameter space in Fig. [I] (namely for f 3 < 0.1 in case of Ms = 10 16 GeV). In these cases, 
the matter-parity symmetry is never restored after inflation, despite its low breaking scale, 
thus avoiding all domain wall problems. 

Secondly, even if the symmetry is restored, domain wall problems can be in principle 
avoided if the discrete symmetry is embedded into an U(l) gauge symmetry, since in that 
case the minima of the potential are continuously connected and the domain walls are 
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Qf Q2 Qt Q% 


-^13 Z14 Z23 Z24 


X 


U(l) fl -x 


11 11 

2 2 2 2 








U(1) (acc.) 
u(1) (acc.) 
u(1) (acc.) 


1 11 1 

2 2 2 2 

1 111 

2 2 2 2 
1111 
2 2 2 2 


-10 1 
0-110 
2 2 2 2 
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Table 4.1: The global symmetries of the superpotential $2.4\) . The standard model fields are 
not charged under U(l)^f c '\ whereas their charges with respect to U(l) B _ L and 11(1)^°°^ 
are the canonical ones. The helds S, S are neutral with respect to the accidental symmetries. 

no longer stable. This mechanism requires however that the gauge symmetry is broken 
after inflation, since only then the domain walls are bounded by cosmic strings and they 
shrink and decay rapidly when the U(l) is spontaneously broken after inflation ends |42] . 
Therefore, in cases where the U(l) B _ L breaking occurs after inflation, our model is free of 
cosmological domain wall issues, despite the low scale breaking. 

The superpotential (12. 4p has four global U(l) symmetries, with charge assignments 
as listed in Table I4TT1 Besides U(1) B _ L it exhibits three additional accidental symmetries, 

U(l)^ acc ' ) , U(l)j ) acc '' ) and the i?-symmetry U(l)^ cc '\ They reflect the absence of linear terms 
and quadratic terms for the fields Zij in the superpotential. These symmetries are either 
explicitly broken by higher order operators, for example trilinear terms which do not in- 
fluence the vacuum ( 12. 5ft . or they are symmetries of the full Lagrangian. In the latter case 
one may gauge the symmetries, since they are anomaly free. However, we do not want to 
have additional massless gauge fields in the model, and hence we rather consider unbroken 
discrete subgroups of these symmetries in that case. 

The accidental i?-symmetry disappears if one introduces mass terms like W = MZ13Z24 
in the superpotential. Such terms are compatible^] with the vacuum under consideration. 
However, if one wants to avoid the prediction of light fermions with masses ~ A 2 /M, 
generated by a see-saw like mechanism in the hidden sector, we need that M C A. This 
may follow if the explicit violation of U(l)^ acc ' ) occurs only via higher- dimensional operators 
in the hidden sector. Alternatively, one may consider that the fields Z^ are moduli of an 
underlying theory, and therefore have no mass terms. Finally, note that the combination 
_R acc - + B — L is orthogonal to the vevs (V12), (V34). After supersymmetry breaking this 
symmetry is broken to its Z 2 subgroup, which is trivial up to hypercharge. 

The absence of linear terms W = C^Z^ is crucial for our model, since such terms 
completely remove the vacuum (12. 5p . It can be explained by promoting a discrete subgroup 
of U(l)^' to a rigorous symmetry of the hidden sector. An important consequence is then 
that the lowest mass states of V13, V14, V23 and V24 become stable together with Z^ due 
to the U(l)^ c ' ) symmetry. As discussed above, this can be problematic since their relic 

8 In addition new supersymmetric minima with (X) = ±A 2 /M occur, where B — L remains unbroken. 
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mass density may overdose the universe, if annihilation into the pseudo NG multiplet is 
not efficient enough. In this cases it would be better to break the U(l) 6 symmetries 
explicitly, to allow the decay of the hidden sector particles. As an example, we note that 
the linear combination U(l) 0+6 alone is needed in order to forbid the linear terms in Z^. 
This combination is preserved if we enlarge (12. 4p to 

W dyQ = X (Pf (V tJ ) - A 2 ) + AV 13 (Z 13 + Z u ) + ...+AV24 (Z 23 + Z 24 ) . (4.1) 

We can then gauge U(l) . b and assume it is spontaneously broken to a Z 2 subgroup at a 
high scale A', such that the heavy hidden sector states have odd parity. If we now identify 
this Wj2 with the Peccei-Quinn like of the SSM (cf. |26| . in SU(5) notation)^ 

P (10) = P(H d ) = - , P (5*) = P (H u ) = P (N c ) = + , (4.2) 

decay channels of the hidden sector states to H u Hd open up. This can resolve the overdos- 
ing issue of the hidden sector from the start. Since the breaking of this is then related 
to the generation of the fi term, we expect that also the linear terms in Zij remain small, 
since they are additionally suppressed by the scale A'. 

5 Conclusions 

In this paper we have discussed a possible mechanism that generates a small matter-parity 
violation in the visible sector. Such a violation of matter-parity can resolve the tension 
between gravitino LSPs and high reheating temperatures as required by thermal leptogene- 
sis. In our scenario, visible sector and part of the hidden sector are simultaneously charged 
under a gauged U(1) B _ L . This XJ(1) B _ L is broken at a high scale Ms by right-handed 
neutrino masses to its matter-parity subgroup. Matter-parity is then subsequently bro- 
ken completely in the hidden sector by a SU(2) hid quark condensate at a scale A. This 
mechanism induces a small vev ~ A 2 /Ms for the right-handed neutrinos, giving rise to 
bilinear matter-parity breaking in the visible sector. Due to matter-parity breaking, the 
NLSP, typically a stau or neutralino, can then decay into standard model particles before 
conflicting with BBN. 

Since the hidden sector is almost separated from the high-energy B — L breaking sector, 
the hidden sector possesses an accidental global B — L symmetry which is explicitly broken 
only through the interactions with the right-handed neutrinos. The effects of the explicit 
breaking are suppressed by I/M5 at low energies, and hence the pseudo NG modes arising 
from the B — L breaking in the hidden sector acquire only very small masses. The tightest 
constraints on our proposed model come from the dynamics of these pseudo NG modes. 
The pseudo NG boson itself acquires O(MeV-GeV) masses, remains practically stable, and 
can be easily overproduced if the coupling between hidden and visible sector is too large. 
The fermionic and bosonic partners of the pseudo NG boson, however, acquire (9(100 GeV) 
masses, decay early in the universe and can reintroduce problems with BBN. We showed 

9 This forbids the dangerous superpotential operators 10 10 10 5*, 10 5* 5*. but not 5*H U . 
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that the model can be phenomenologically viable for large enough B — L breaking scale 
Ms > (9(10 13 GeV), and for condensation scales A ~ 0(1O 3 -1O 8 GeV), when the coupling 
to the lighter right-handed neutrinos is suppressed. 

An important consequence of gravitino dark matter scenarios with small breaking of 
matter-parity is that the gravitino decays, producing gamma-ray lines which are potentially 
observable with instruments like Fermi LAT. We found that within our model, large values 
of the matter-parity violation are favored, since this reduces the lifetime of the problematic 
pseudo NG modes. This suggests lifetimes of the gravitino close to the experimental limit, 
as well as short decay lengths of the NLSP. For a stau NLSP, we found a lower limit for 
typical stau decay lengths that is given by > (9(100 m). Such events are in principle 
detectable at the LHC 
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A Fermion Mixing Matrix 

The superpotential contributions (12.21) . (12.41) and (12. 6p imply mixings between the 13 neu- 
tral fermions 

B,W,hlh° d ,u h N^V 12 ,V u ,X , (A.l) 

where the first four fields are the MSSM gauginos and neutral higgsinos, z/j and Nf are 
the left- and right-handed neutrinos, respectively, and the last three fields are fermions of 
hidden sector chiral multiplets. In this basis the mass matrix is to leading order (and up 
to order one coefficients) 
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where £ = A/Mp and v u denotes the vev of the bosonic component of the superfields H u , 
and we took / 3 = 1 and f 1 — f 2 — for definiteness. The coefficients ctj and are non-zero 
for non-vanishing vevs of the the left- and right-handed sneutrinos, respectively. 

The 4 x 4, 6 x 6 and 3x3 blocks on its diagonal correspond to the MSSM neutralinos, 
the left- and right-handed neutrinos and three hidden sector fermions. In the basis in 
which the neutralino and the hidden sector mass matrices are diagonal and the left and 
right-handed neutrinos decouple, the matrix reads (only leading order terms are shown) 
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a ai A 2 /M s 


/3 ai A 2 /M s 
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(A.3) 

with neutrino mass matrices M Vhij = -vlJ2k^ik^kj T '/ M k, M URij = Mrfij + 0(vl/M s ), 
and diagonal neutralino masses M xa/3 = m Xa 5 a/ 3, a, (3 = 1, ... ,4. 

We deduce that the left-handed neutrinos mix with lightest neutralino Xi an d the light 
pseudo NG fermion ^ as 

L ' y/2M s m^~ ' ^10^^ VlO 6 GeVy 1 VlO 16 GeVy 1 V100 GeV7 ' 

(A.5) 

B Pseudo NG mode masses 

We will shortly discuss how the pseudo NG modes of B — L breaking in the hidden sector 
acquire masses after SUSY breaking. To this end we will assume that SUSY breaking is 
only mediated by SUGRA effects to the matter-parity breaking sector, and that the Kahler 
potential is canonical. 

The mass-generation of the fermionic and bosonic partner of the pseudo NG boson, ip 
and a, is due to SUSY breaking and can be understood when considering the hidden sector 
alone, neglecting couplings to the visible sector. After breaking of SUSY, soft terms are 
induced in the scalar potential of the hidden sector. They have the form 

Koft = ml/2 (|AT + |U 12 | 2 + |U 3 4| 2 ) + 2m 3/2 (A 2 X + h.c.) , (B.l) 
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where we suppressed 0(1) prefactorsl^l The supersymmetric part of the scalar potential 
reads 

VW = \V 12 V M - A 2 | 2 + \XV l2 \ 2 + \XV M \ 2 , (B.2) 



and (IB.ip and (IB.2|) have together a minimum ax^_ 



1 + 2^ J . (V a ) = (V U ) * A - -j^ • (B.3) 

The vev of X induces a Dirac-mass term for the fermion partners of V\ 2 and V34, which 
finally gives a mass close to m 3 / 2 to the fermionic partner of the pseudo NG boson, if). 
By expanding the fields like V\ 2 = (Vi 2 ) + p/2 + . . . and V34 = (V34) — p/2 + . . . one can 
furthermore derive that the mass of the bosonic partner of the pseudo NG boson is given 
by m p ~ 4m3/2- For definiteness, we use these masses in the phenomenological analysis, 
and assume that ip is somewhat heavier than the gravitino 1/^3/2, although the exact values 
can change depending on the details of the Kahler potential and SUSY breaking sector. 

Without coupling to the visible sector, the pseudo NG boson in the hidden sector 
would remain massless after SUSY breaking. However, the coupling to the right-handed 
neutrinos generates an additional term (/| A 2 /M 3 ) V34V34 in the superpotential, which after 
SUSY breaking give rise to terms like (f 2 rri3/ 2 A 2 /M3) V34V34 + h.c. in the scalar potential. 
Noting that V 34 = A e ~ ia / 2A + . . . , it follows that the mass of the pseudo NG boson in the 
hidden sector is of the order of m a ~ ^3 ^/m 3 / 2 A 2 /M3, see Eq. (12 .9p . 
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